OFF-LINE FEED RATE SCHEDULING FOR REDUCTION OF MACHINING TIME 
AND ENHANCEMENT OF MACHINING ACCURACY IN CNC MACHINING 
BACKGROUND OF THE INVENTION 

(a) Field of the Invention 

The present invention relates to CNC (computer numerically controlled) 
machining technology, and more particularly, to off-line feed rate scheduling for 
reduction of machining time and enhancement of machining accuracy in CNC 
machining. 

(b) Description of the Related Art 

Research into metal machining technologies has taken on increasing 
importance as design and production cycles of a product become shorter in response 
to changing consumer demand. This is particularly true in the present-day industrial 
situation where a company's survival is closely linked to its ability to remain 
economically competitive. Modern production workers perform machining by relying 
upon empirical machining conditions even after a tool path is generated using CAM 
(computer aided manufacturing) software. Production workers depend mainly on stable 
machining conditions such that production time is reduced. Further, in transient cut 
regions such as a corner section where a tool path is encountered, machined surface 
errors occur severely because production workers randomly adjust a feed rate in these 
regions. In addition, many processes are required to remove these machined surface 
errors generated during NC machining. 

These problems may be overcome through the systematic selection of 
machining conditions. Since a feed rate scheduling system gives a manufacturing 
machine user optimum machining conditions that take into consideration ease of 
machining, machining stability, and machined surface errors, it can be a useful tool in 
many respects (e.g., to reduce production costs and to make more various kinds of 
small production). 

In CNC machining, although the need to optimize feed rate, which is related to 
productivity improvement, is becoming increasingly important, it has not yet become a 
common technology that can be used for the manufacture of a wide variety of products. 
This is a result of the fact that analysis is not definitely performed with respect to the 
cutting process such that it is difficult to anticipate machining results of various cutting 
conditions. In order to determine optimum cutting conditions for CNC machining, a 
process of trial and error is undertaken. Since it is difficult to format this knowledge for 
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input to a database, a new production worker must start anew until he or she becomes 
familiar with the processes. 

CAM software companies are introducing feed rate optimization methods to 
improve their competitiveness. However, feed rate optimization methods used in most 
5 CAM software available today uniformly maintain an MRR (material removal rate) in 
consideration of a geometric relation between a workpiece and a tool. However, in 
methods utilizing MRR, it is not possible to put precision adjustment of a cutting force 
and a cutting force coefficient in database form. Since the intensity of the cutting force 
may vary greatly according to changes in spindle speed even with the same MRR, 
10 such methods result in reducing competitiveness. Accordingly, a feed rate scheduling 
method that is based on a cutting force model is needed. If accurate, the cutting force 
model would be able to be effectively applied to feed rate scheduling. 

SUMMARY OF THE INVENTION 
In one exemplary embodiment of the present invention, there is provided an 
15 off-line feed rate scheduling method that improves machining precision and productivity 
during CNC machining by developing an off-line feed rate scheduling model that 
optimizes feed rate in CNC machining based on a cutting force model and a surface 
error model. 

In an exemplary embodiment of the present invention, an off-line feed rate 
20 scheduling method of a CNC machining process includes selecting a constraint 
variable and inputting a reference value related to the constraint variable; estimating a 
cutting configuration where a maximum constraint variable value (CW) occurs through 
ME Z-map modeling; receiving the estimated cutting configuration and estimating a 
specific rotation angle (<}) s ) where the maximum constraint variable value occurs 
25 through constraint variable modeling; calculating a feed rate that satisfies the reference 
value of the constraint variable at the estimated specific rotation angle; and applying 
the calculated feed rate to the NC code. 

The cutting force or machined surface error may be selected as the constraint 
variable depending on machining conditions. 
30 Calculating the feed rate includes inputting specific feed rates U and f 2 (fi < f 2 ); 

calculating maximum constraint variable values CW, and CW 2 corresponding to the 
feed rates f1 and f2, respectively, at the specific rotation angle; approximating a feed 
rate f nex t that corresponds to a reference value RV of a constraint variable value using 
the formula 
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f f , (m T -cw x yj 2 - /,) 

calculating a constraint variable CW next in the case where the feed rate is f ne xt; 
and determining using the formula below if the constraint variable value CW ne xt when 
compared to the reference value RV is less than an error limit, applying the feed rate 
5 fnext to the NC code when it is less than the error limit, replacing the feed rate f 2 by f nex t 
and repeating the process of obtaining f ne xt when this value is not less than the error 
limit and the reference value RV is greater than the constraint variable value C\A/ ne xt, 
and replacing the feed rate fi by f nex t and repeating the process of obtaining f ne xt when 
this value is not less than the error limit and the reference value is not greater than the 
10 constraint variable value CW next , 

where 

CW next -RV 

— < Error Limit 

RV 

The computing a cutting configuration through ME Z-map modeling includes 
searching for node points located in cutting area; identifying whether a target node is 
15 the edge node (the node closest to cutter edge) or not; calculating and updating a 
height value of each node in cutting area; moving a target node if it is an edge node 
and storing movement direction angles; and computing the cutting configurations using 
the stored angles. 

If a difference between a distance from a tool center to a target node and a 
20 tool radius is smaller than a movement limit, this node is designated as an edge node. 

Also, the cutting configurations computed through ME Z-map modeling include an entry 
angle, an exit angle, an axial depth of cut, and so on. 

BRI^F DESCRIPTION OF THE- DRAWINGS 
The accompanying drawings, which together with the specification, illustrate 
25 exemplary embodiments of the present invention, and, together with the description, 

serve to explain the principles of the present invention. 

FIG. 1 is a schematic diagram showing a structure of an off-line feed rate 
scheduling system according to an exemplary embodiment of the present invention. 

FIG. 2 is a schematic diagram showing the formation of a flat end milling 

30 process. 

FIG. 3 is a drawing used to show a tool rotational angle and an angle at a 
position to be cut at a cross section of a flat end milling tool. 
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FIG. 4 is a drawing used to describe a method of calculating a cutting force 
coefficient. 

FIG. 5 is a flow chart of a method for calculating a cutting force coefficient. 
FIG. 6 is a drawing used to describe a method of calculating a thickness of an 
uncut chip. 

FIGS. 7A and 7B are schematic views showing a shape and coordinates of a 
ball end milling tool. 

FIG. 8 is a drawing used to describe a size effect according to the geometry of 

a tool. 

FIG. 9 is a schematic view showing a local helical angle 9 h! and a local 
effective helical angle 6 h i c in a ball end milling tool. 

FIG. 1 0A is a schematic view showing feed rate components in 3D machining 
of ball end milling, and FIG. 1 0B is a drawing used to describe the calculation of an 
uncut chip thickness in ball end milling. 

FIG. 11A is a drawing used to describe a first case in which machined surface 
errors are generated, FIG. 11B is a drawing used to describe a second case in which 
machined surface errors are generated, and FIG. 11c is a drawing used to describe a 
third case in which machined surface errors are generated. 

FIG. 12 is a drawing used to describe a machined surface error. 

FIG. 13 is a drawing used to describe run out of a tool. 

FIG. 14 is a drawing used to describe tool deflection. 

FIGS. 15A and 15B are schematic views showing a shape of cusps, and FIG. 
15C is a drawing used to describe a method of calculating a cusp height. 

FIG. 16 is a drawing used to describe a fundamental notion for ME Z-map 

FIG. 17A is a drawing showing a edge area including a edge node, and FIG. 
17B is a drawing used to describe the moved edge node. 

FIG. 18 is a drawing used to describe the movement limit of an edge node. 

FIG. 19 is a schematic view showing the case where there are two or more of 
each entry and exit angles. 

FIG. 20 is a drawing used to describe grouping method of disk elements. 

FIG. 21 is a flow chart of an algorithm for cutting configurations calculation. 

FIGS. 22A and 22B are drawings used to describe the determination of a 
reference cutting force based on a transverse rupture strength of a tool. 

FIG. 23 is a flow chart of an algorithm used to calculate a feed rate that 



satisfies a reference value in order to adjust a constraint variable in an NC code 
interval. 

FIG. 24 is a graph showing acceleration and deceleration characteristics of a 
machine tool. 

5 DETAILED DESCRIPTION 

Exemplary embodiments of the present invention will now be described in 
detail with reference to the accompanying drawings. 

FIG. 1 is a schematic diagram showing a structure of an off-line feed rate 
scheduling system according to an exemplary embodiment of the present invention. 

10 An off-line feed rate scheduling system includes a cutting force model, a 

machined surface error model, an ME (moving edge node) Z-map model, and an 
optimum feed rate scheduling algorithm. 

As shown in FIG. 1, after selection of an NC code, machining information, and 
variables for feed rate scheduling, that is, constraint variables (cutting force and 

15 machined surface errors), feed rate scheduling is performed when a reference value, a 

simulation interval, and an acceleration/deceleration time are input. In the case where 
the cutting force is established as the constraint variable, the reference value is 
automatically calculated without requiring additional input; and in the case where the 
machined surface error is established as the constraint variable, the user inputs a 

20 desired machined permissible error as the reference value. 

Further, cutting configurations such as a tool traveling direction, a tool 
entry/exit angle, an axial depth of cut, etc. are calculated from the ME Z-map and 
transmitted to the cutting force model and the machined surface error model. Next, 
cutting force, tool deflection, and machined surface errors are calculated from the 

25 cutting force model and the machined surface error model. So that these values satisfy 
reference values, the feed rate scheduling algorithm is used to find the desired feed 
rate and provide an optimized feed rate to the NC code. 

The method of performing feed rate scheduling varies depending on the 
constraint variable. In this exemplary embodiment, cutting force or machined surface 

30 errors may be selected as the constraint variable. That is, in the case where the user 
intends to adjust cutting force or machined surface errors, a feed rate scheduling 
method according to the selection of each constraint variable may be used. The user 
may select cutting force or machined surface error as the constraint variable according 
to the step of rough cutting, semi-finish cutting, and finish cutting. Typically, cutting 
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force is selected as the constraint variable in rough cutting and semi-finish cutting, and 
machined surface error is selected as the constraint variable in finish cutting. 

In the following, a method of calculating and obtaining models of each of the 
cutting force and machined surface error, and of calculating a cutting configuration 
5 through the ME Z-map model, is presented. Also, a method is described in which these 
constraint variable (cutting force and machined surface error) models and the ME Z- 
map model are used to calculate an optimum feed rate according to allowed reference 
values. 

In the case where cutting force is selected as the constraint variable, a cutting 
10 force model is required to estimate the cutting force. In this exemplary embodiment, 

cutting force models with respect to flat end milling and ball end milling have been 
calculated. 

[Flat end milling cutting force model] 

FIG. 2 is a schematic diagram showing the geometry and the coordinate 
15 system of a flat end milling process, and FIG. 3 is a drawing used to show cutter 

rotation and cutting edge location angles on a cross section of an end-milling cutter. 

The end milling cutter, as usual, is divided into a finite number of disk 
elements and the total x-, y-, and z-force components acting on a flute at a particular 
instant are obtained by summing up the force components acting on each individual 
20 disk element. 

In FIG. 3, each angle position of an i-th tooth and a k-th axial direction disk 
element in an initial state (when tool rotation is not considered) is obtained using 
Formula 1. 

[Formula 1] 

a(i, k) = (i— l)<f> c + (kAa + /la/2) 

25 

where i is a cutter tooth index, k is an index of a z-axis disk element, z7a is a 
height of the z-axis disk element, 0 h is a helix angle, and R is a radius of the cutter. 

In Formula 1, an angle position of a cutter edge is obtained using Formula 2. 
[Formula 2] 

<f>(zj,k) =a(z 9 k) + 0U) 

eu) =-jde 

30 

where j is an index of a cutter rotation angle, 6 (j) is a cutter rotation angle, 
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and /I 6 is the amount of increase in the rotation angle. 

If the cutting force components on rake surface is divided into a normal 
pressure force dF n ( 4> ) and a frictional force dFf( $ ), Formulas 3 and 4 are obtained. 

[Formula 3] 

dF n (a,0) = K n A(0)n(a)dA c 

[Formula 4] 

dF/ia,0) = K / K n A(0)~T\(a)dA c 
where 

dA= t c {cf>) cos a r {-^-) 

cos U " , and 



A{6) = 



cos 6 — sin 6 0 
sin 6 cos 6 0 
0 0 1 



In the above conditions, a r is a rake angle and t c (<t>) is uncut chip thickness, 
which are calculated considering tool deflection and run out. 

If the cutting force is shown in a rectangular coordinate system using the 
normal pressure force and the frictional force, Formula 5 is obtained. 
15 [Formula 5] 

F x ( i, j, k) = [C i K n cos(<f> - a r ) + K f K n C 3 cos <f> - K f K n C 4 sin ( <f> - a r )]t c ( <f>) B { 

F y (iJ, k) = [ C,K n sin(^- a r ) + K f K n C,sm<f> + KjK n C 4 cos(<f>- a r )]t c (<f>)B x 

F z (zJ,k) =[-QK n + K / K n C 3 ]t c (frB l 

where C1 and C2 are constants, and are calculated by the following: 

__ cosd h _ sm6 k 

Li — : — a — C 2 - — : — n — * COSa r 

and C3, C4, and C5 are values determined by chip flow angle 9 C and is 
20 calculated as in the following: 

C 3 = sin 0 h { sin 0 C — cos 0 c cot 0 tk ) 

cos0 r 



Ci ~ sin6 tk 
C 5 = cos 6 h ( sin 0 c — cos 0 c cot 6 tk ) 



l 



where cas0 th = sina r - sind k 

Accordingly, cutting force components of each direction according to the 
rotation angle of the cutter are obtained by Formula 6. K nt K f , and 0 C according to 
cutter rotation angle can be obtained using F Xj F y , and F z measured through 
experimentation. 

[Formula 6] 



FIG. 4 is a drawing used to describe a method of calculating a cutting force 
coefficient, and FIG. 5 is a flow chart of a method for calculating a cutting force 
coefficient. 

The cutting force coefficient may be obtained using an instantaneous cutting 
force generated by a cutting edge of a tool disk element obtained by the method as 
shown in FIG. 4. Flat end milling includes a bottom edge and a helical edge, and since 
the cutting force model is realized with respect to the geometry of the helical edge, the 
cutting force coefficient is obtained using the cutting force generated by the helical 
edge. 

With reference to FIG. 5, if a cutting force measured under conditions of a 
small axial depth of cut is subtracted from a cutting force measured under conditions of 
a large axial depth of cut, a cutting force generated by a helical edge of a unit disk 
element can be obtained. That is, cutting forces F 2 and Ft are measured respectively in 
axial depths of cut d a2 and d a i (d 3 2 > d a i). Next, if F 1 is subtracted from F 2 after 
synchronization of each cutting force measured value and estimated value, cutting 
force F(ij,k) applied to the helical edge of a unit disk element can be obtained. Using 
this cutting force F(ij,k), tool runout and cutting force coefficient can be obtained using 
formulae 7, 8, and 9. The developed cutting force coefficient appears in the following 
independent form with respect to cutting conditions. The non-linearity is a result of the 
fact that the cutting force coefficient is influenced by size effect in an area where uncut 
chip thickness reduces. 

[Formula 7] 





FJU) 
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\n(K„(i,j,k)) = A -(4 - A 2 )e-^'- J - t>rU 

[Formula 8] 

K f (i,j,k) = B X - (B, - B 2 )e-^> J ' k))BA 

[Formula 9] 

e c (i,j,k) = C l -(C l -C 2 )e- (C "^ t » c ' 

5 

FIG. 6 is a drawing used to describe a method of calculating uncut chip 
thickness. 

An uncut chip thickness model is needed to calculate cutting force in Formula 

5. 

10 Formula 10 is used to model uncut chip thickness using the cutter geometry. 

This model has the advantage of requiring only a minimal time to calculate uncut chip 
thickness, but it is limited in its precision as it simply reflects a geometric state. 
[Formula 10] 

t c (<t>) = f t sxncf> 

15 In this exemplary embodiment, uncut chip thickness is calculated while 

following the movement of the tool. That is, a difference between a track of a previous 
edge left on a surface of the workpiece and a present position of the tool is used to 
obtain the uncut chip thickness. As shown in FIG. 6, in order to obtain the uncut chip 
thickness, it is necessary to first know x am and y am - x tm anc ' ytm are determined using 

20 Formulae 1 1 below. x am and y am in the Formulae 1 1 are a center position of the tool 

when a previous edge (an edge in m=1~n f ) passes an angle position (4>) of a present 
edge, and x tm and y^ indicate positions of a previous edge at x am and y am . 
[Formula 11] 

y, m =(-l/tan <f>) • x tm 

(x tm -x am ) 2 + (y tm -y a J 2 =R 2 

25 From Formula 11, the uncut chip thickness may be calculated using Formula 

12. 

[Formula 12] 
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tcmiij, k) = R-\J x\ m + y l tm 
t c = Max [ 0 , Min ( t cm ( i 9 /,*))] 

In the following, cutting force model, cutting force coefficients, and uncut chip 
thickness model in the case where a ball end milling tool is used are presented. 
[Ball end milling cutting force model] 

FIGS. 7 A and 7B are schematic views showing a geometry and coordinates of 
a ball end milling tool. 

The envelope of the ball part of the ball-end mill can be expressed by 
Formulae 13 and 14 below: 

[Formula 13] 

[Formula 14] 

R(zf =x 2 +y 2 

where R Q is a ball radius at a circle of a hemisphere center, and R(z) is a 
radius in an x-y plane at a specific z. 

A local helical angle is determined by Formula 15 below. 
[Formula 15] 

<!.(»). t M -( * w B ' m *» ) 

In FIG. 7A, a vertical unit vector (^) and a chip flow vector (%) are defined on 
the rake surface of a differentiated cutting edge. As a result, the normal pressure force 
and the frictional force acting on the rake surface can be derived. In order to derive the 
formula for ball end milling, the unit vectors are obtained when the cutter rotation angle 
0=0 and the local helical angle 9 h =0. The unit vectors having rake angle (a r ) are 

designated as follows. That is, the unit vectors n normal to the rake surface are 
determined by Formula 16 below. 
[Formula 16] 

« = cosa r / -sxnaj 

The unit vectors tangent to the cutting edge are obtained by Formula 17 below. 
[Formula 17] 
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-/ x dyj + dzk - /, -t 1 

jdy 2 +dz 2 f 2 A 



dr 



k 



j _ dy_ _ (z— R 0 ) 

dz ^R 0 2 - (z-R 0 ) 2 f 2 = y/ff+l* 

where and 

The unit vector b that is normal to the tangential direction of the cutting edge 
and on the rake surface is as shown in Formula 18 below. 
5 [Formula 18] 

r - - 1 . - 1 -r /, 

b = t x n - — sin a i h cosa / - — cosa r A: 

h h h 

Accordingly, the chip flow vector is obtained by Formula 19 below. 
[Formula 19] 

f c = cos0 c b(a) + sin0 c f(a) 

= cos0 c (— sin a r i + —cosaJ-—cosa r k) + sin0 c (—j + —k) 

J 2 J 2 J2 J2 J 2 

- cos0 c sin a r —i + (cos0 c — cosa r + sin 0 C — ) j + (-cos# c — cosor r + sin 0 c —)k 

fi fi y^> fi 

10 The cutting force on the rake surface is determined by two orthogonal 

components, that is, a normal pressure force and a frictional force as shown in 
Formulae 20 and 21 below. 
[Formula 20] 

dF n (iJ,k) = K n T(^0 M )ndA c 

15 [Formula 21] 

dF f (iJ,k) = K f K n T(0,0 M )f c dA c 

The transformation matrix resulting from the cutter rotation angle 9 and the 
helical angle 9 h is obtained using Formula 22 below. 
[Formula 22] 



cos(-0 + a) -sin(-0 + a) 0 

sin(-0 + a) cos(-0 + a) 0 
0 0 1 

cos (j) cosOfj -sin^ cos </> sin 6 h 

sin <p cos 0 h cos <t> sin <f> sin ^ 



cos(0„) 
0 

-sin(^) 



0 sin(0») 

1 0 

0 cos(0„) 



- sin 0 h 



COS0,, 



If the normal pressure force and the frictional force on the rake surface are 
converted into a cutting force of three axial directions, the following Formulae 23, 24, 
and 25 result. 

[Formula 23] 



FJiJ,k) = [A^„(cosa r cos <t> cos 0 h + sina r sin^) + A^ n AT / [cos^ c sin or r — cos ^ cos ^ 

1 f f \ 

- (cos# c — cosa r + sin 0 C — )sin <j> + (-cos# c — cosa r + sin 0 C — )cos^sin0 A ]]d£4 e 

J 2 f2 f2 fi 

= [AT fl (cosa r cos <f> cos 0 h + sina r sin <p) 

1 l / 
+ K n K f cos 0 C [sin a r — cos^cos0 A cosa r sin^ -cosa r cos ^ sin 0 h ] 

J 2 J 2 J2 



+ K n K f sin # e [- — sin ^ + — cos 0sin 0 h ]dA c 



[Formula 24] 



F v 0J>k) = [K„( cosa r sin ^ cos 0 ^ -sin a, cos^) + A^A^cos^ sina r — sin <j> cos 0 h 

J 2 

\f f I 

+ (cos 0 C — cos a r + sin # c — ) cos ^ + (— cos 0 C — cos a r + sin 6? c — ) sin <j> sin 0 A 

y*2 fi /i y*2 

= [A^(cosa r sin <ficos0 h -sin a r cos^) + 

+ K n K f cos# c (sina r — sin^cos# A + — cosa r cos</>-— cosa r sin^sin# A ) 

/2 -Z2 /2 

/* 1 

4- K n K f sin# c (— cos^ + — sin^sin^ )]ci4 c 

/2 ^2 



[Formula 25] 
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F : (iJ,k) - [-K n sin0 h cosa r + K n K f [- cos 6 C sina r — sin^ 

fi 

+ (-cos# c — cosa r + sin# c — ) cos 0 h ]cIA c 
fi fi 

1 / 

= [-K n sin0 h cosa r + K n K f cos# c (-sina r — sin^ A cosa r cos0,,) 
+ K n K f sin0 c -^cos0 h ]]dA c 

J 2 



Each directional cutting force component according to the rotation angle of the 
cutter can be obtained using Formula 26 below. Using the cutting force Fj measured 
from experiment, the cutting force coefficient according to cutter rotation angle, K n , K f , 
and 0 C can be obtained. 

[Formula 26] 



ijk 



If Formula 26 is put in matrix form, Formula 27 is obtained. 
[Formula 27] 



F x 








F. 





A n A ] 



*21 



A 



31 



12 
^32 



A 



13 



*23 



A 



33. 



K l 



K 



3 J 



where 
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K 2 =cos0 e K n K / 



K 3 =sm0 c K a K f 



Ai = B\ y^y^ ( coscg r cos^cos# A + sin<z r sin^) ■ f c (0) 
a i 

1 1 f 

A2 = AZZ( s * na ' — cos ^ cos ^ A -— cosa r sin^ - ^f-cosa r cos^sin^) t c ($) 

k 1 72 /2 

^13 = B iI!Z(' : 7 sin ^ + y C0S ^ sin ^)'^^ 

* r 

11/" 

^22 ~ ^ y^y^ C 5 " 1 ^ — sin 0 cos 0^ h cosa r cos^ L cosa r sin^sin^) -t c (<p) 

k i f 2 fz ft 

f 1 

^23 =5 iXZ(y cos ^y sm ^ sm ^)-^W 

k i Ji J2 

' 4 3i = 5iZI(-sin^cosa r )'r c W 

* 1 

^32 = Z Z ( _sin °> ~r sin0 h- -r cos a r cos # A ) • l c 

K 1t K 2l and K 3 can be obtained from the measured cutting force, and the 
cutting force coefficient can be calculated using Formula 28 as follows. 
[Formula 28] 

0 c= tan-'A) 

_ K 2 
f " cos^ 



FIG. 8 is a drawing used to describe a size effect according to geometry of a 



tool. 



FIG. 8 shows the changes in the length of the cutter edge per unit disk 
element. If the height of a disk element (Aa) is constant along the z-axis, the length of 
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the cutter edge (ds) increases rapidly as the z position of the disk element approaches 
zero. This change in the cutter edge length alters the area of the flank face contacting 
the workpiece accordingly. The force increases as the z position approaches the 
bottom disk, although the uncut chip thickness remains unchanged at the same cutter 
edge location angle. This directly influences the size effect of a ball-end mill in the form 
of the plowing force acting on the rounded edge or flank face of the cutter. In the 
present invention, the disk edge length is used to obtain this effect such that the cutting 
force coefficient changes in the axial direction are considered. 

In ball end milling, after obtaining a cutting force generated by disk elements 
by the method shown in FIG. 8, the cutting force coefficient is obtained. The cutting 
force coefficient is determined using Formulae 29, 30, and 31 below. Since the cutting 
force coefficient formulae of the present invention are independent of cutting conditions, 
cutting force coefficients can be estimated through two calibration tests. As shown in 
Formula 29 below, K n is obtained considering size effect generated by each disk edge 
and with respect to rescaled instantaneous uncut chip thickness t CT using each disk 
edge length d s . 

[Formula 29] 



\n(K n ) = A l -(A l -A 2 )e-^ )44 



t c (i,j,k)/ co$0 M 



ds{k)l Aa 
where 

Since K f and 0 C are not affected by the length of each disk edge, they are 
obtained using the instantaneous uncut chip thickness t cr . 
[Formula 30] 



K f =B x -(B,-B 2 )e-^ 



[Formula 31] 



C -C 
= ^! 1_2 + 



° c ° hlc ~iHt c /c 3 f* +C2 



where 0 h i c is a local effective helical angle of a ball region that affects chip 
flow. The localized helical angle 9 h i and the localized effective helical angle 9 h i c 
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appear in FIG. 9. 

FIG. 1 0A is a schematic view showing feed rate components in 3D machining 
of ball end milling, and FIG. 1 0B is a drawing used to describe the calculation of uncut 
chip thickness in ball end milling. 

In the 3D machining shown in FIG. 10A, the uncut chip thickness is calculated 
using Formula 32 below. 

[Formula 32] 

t c (z)= f t • sin <fi • cosy + DR(z) 

where DR{z) = R{z) - R{z + f r sin y/) and thjs term refers to the 

generation of an uncut chip thickness by changes in a tool disk radius at the same 
height by axial direction feed component in three-dimensional machining. In Formula 

32 if f c( z ) =R( < z ) z + f t siny <0 

In the 3D machining shown in FIG. 10A, in order to take into consideration 
changes in uncut chip thickness caused by vertical feed component, changes in radius 
according to cutter axial direction as shown in 33 below were considered. 

[Formula 33] 

-1 

tan^ 

(Xfin ~ X am + X a f + <X, " y am + ^ f = R<J 

where Ram = R(z + mf t sin$) 

Accordingly, uncut chip thickness is calculated using the following Formula 34 
in Fig. 10B. 

[Formula 34] 

tern ('"> h k > ™) = R ( Z ) ~ 4 X imi m f + ynni m f 

t c (z) = Max[Q, Minify (/, j, k, w))] 

where tcmflj.k.nn) refers to a possible chip thickness, x tm and y tm indicate a 
cutting edge position (m=1~n f ) at a previous edge, and t c is an actual uncut chip 
thickness. 
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[Machined Surface Error Model] 

In the following, a machined surface error model is described. This model is 
used to estimate surface errors when machined surface errors are selected as the 
constraint variable. 

As shown in FIG. 11 A, if it is assumed that ideal movement occurs only by a 
nominal feed rate, a trace of a disk element can be assumed to be a circular arc. Such 
an assumption is suitable when the cutting feed rate is extremely small compared to 
the diameter of the tool. It is used by many research scientists. In the present invention, 
the machined surface error model for considering tool runout and deflection is 
developed from the cutting force model. During machining, even if a center coordinate 
of a disk element strays from a nominal value by tool deflection and runout, the 
movement as shown in FIG. 11B occurs. Further, in FIG. 11c, there are shown 
machined surface errors generated by ball-end mill geometry and a tool path distance. 
In the present invention, machined surface errors occurring in the cutting process are 
calculated by estimating these machined situations. 

Surface error estimation in first case (FIG. 1 1 A) 

Cusp errors (or cutter marks) that may develop by flat end milling and ball end 
milling are generated by feed per tooth f t and tool radius R as shown in FIG. 12. The 
cusp error C h of FIG. 1 1 A is calculated using Formula 35 below. 

[Formula 35] 



Surface error estimation in second case (FIG. 11B) 

Tool runout generated during tool set up may be determined by offset and its 
angle term as shown in FIG. 13. Variations in tool center coordinates by tool runout are 
obtained by Formula 36 below. 

[Formula 36] 



C k =R-A 





M 



x p (j) =—p sin (a ntt- 60*)) 
y P (j) =- pcos(a nm - d(j) ) 

FIG. 14 shows tool deflection when it is assumed that the tool is a two stepped 
cantilever beam and a static force is applied to the tool. A formula used to obtain a 
position of each disk element using the force applied to each z-axis disk element is 
Formula 37 below. 

[Formula 37] 

[<L f -L p -z> 3 -{L f -zf + 3-L, -(L, -z) 2 ] 



6EI f 



6EI f 



[ <Lf -L p -z>* -(L f -zf+3.L p -{L f -zf) 



where L s and L f are the shank part length and the flute length, respectively. 
L p indicates the length from the position where the cutting forces act to the boundary 
10 between flute and shank part. E is Young's modulus. The moments of inertia for the 

shank and that of the flute part were calculated as 7 / =Z) 4 /52.84 and I s =7r -D 4 1 64, 
respectively 

Using the cutting edge coordinates obtained while determining uncut chip 
thickness according to tool movement, machined surface errors are more easily 

15 estimated. In the case of actual machining, machined surface errors may be obtained 
using cutting edge position (p xk) p yk )i at the instant when the cutting edge position 
angle, which passes the machined surface, is 0. Further, using (p xk , p y k )i coordinates at 
each disk element, three-dimensional machined errors may be formed. The (p Xfk , p y , k )j 
coordinates show the passed positions of the actual tool edge, that is, the positions 

20 including machined surface errors. These machined surface errors may be calculated 
using the following Formula 38. 
[Formula 38] 

E(j\k) = y p (j)+Y d U f k) 
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Surface error estimation in third case (FIG. 11c) 

FIGS. 15A and 15B are schematic views showing a shape of cusps, and FIG. 
15C is a drawing used to describe a method of calculating a cusp height. In FIG. 14, a 
distance between tool centers varies by a tool path distance TPD and edge feed rate f t , 
and this affects cusp height. The greatest cusp height is formed when the distance 
between tool centers is at a maximum. Cusp height is calculated by Formula 39 below. 

[Formula 39] 

D = yj(TPD 2 +f t 2 ) 

c fc = *-^* a -(f) 2 

[ME Z-Map Model] 

As illustrated in FIG. 16, a fundamental notion for ME (Moving Edge node) Z- 
map is to move the edge node toward the boundary of the cutter movement. The edge 
node moves in the direction of 6 m which is the angle of the edge node from the cutter 
center, so that the coordinate of the node (i, j) is changed to: 

[Formula 40] 

x = i g + a x 

y = j-g + <*> 

where g is the elemental grid size, and a x and a y are the distances moved from 
the original coordinates shown in FIG. 16. The value of a x and a y should be smaller 
than the elemental grid size. 

In the simulation process, every node position is calculated one after another 
and the node located in edge area is defined as the edge node as illustrated in FIG. 
17A. FIG. 17B shows that the boundary of the swept volume can be described 
precisely by moving the edge nodes are designated only in the boundary region of the 
swept volume; while the remaining nodes are used to describe the other regions. 

The procedure of moving an edge node consists of the next two steps. The 
first step is to identify whether a target node is located in the edge area or not. If the 
difference between the cutter radius and the distance between the target node and 
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cutter center is smaller than the movement limit of an edge node as illustrated in FIG. 
18, then the target node becomes the edge node. The movement limit can be 
represented by Formula 41. 
[Formula 41] 



KJ n = 



-a x cos0 m ±g|cos0 m 



cos 2 6 m 
-a y sin G m ±g|sin 0 m \ 



for C > 0 
for C<0 



2 n ■ 

sin 0 m 

V m 

where 

C= \a x +l o cos0 m \-\a y +'oSin0 M | l o = g/max(| cos0 m |,| sin0 m |) . 

. I 0 is the 

movement limit of the initial node position (when a x and a y are 0), and l p and l n are the 
movement limit for positive and negative directions, respectively. The second step is 
10 moving the target node if it is located in edge area. The moving of node can be realized 

by adjustment of a x and a y values. Then, the updated values, a x ' and a y \ can be 
described in Formula 42. 
[Formula 42] 

<=<* x +(r<-d m yco&e m 

< = a,+(r<-d m )-sm0 m 

15 where r d is the boundary radius and d n is a distance between target node and 

cutter center. 

In order to estimate the cutting force of the milling process, it is necessary to 
calculate a cutting configuration through simulation. The cutting configuration generally 
includes entry angle, exit angle, and axial depth of cut. The computing process of 
20 entry/exit angle for each node is as follows: 

First, it is determined if a target node is located in the cutting area. Second, if 
the target node is the edge node, it is moved and a height value of node is changed. 
Third, if the target node is the edge node, an angle of movement direction is stored. 
Fourth, the entry angle and the exit angle are calculated using the stored angle of 
25 movement. 

In a cutting condition where only one of each of an entry angle and an exit 
angle are existed, the largest value of the stored movement direction angles is set as 
the entry angle, and the smallest value is decided on the exit angle. However, as 
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illustrated in FIG. 19, if there are two or more entry angles and exit angles, 0 n -i is set 
as the exit angle and 9 n is decide on the next entry angle when the interval between 
stored movement direction angles satisfies the condition of Formula 43. 
[Formula 43] 



9. -0. . >2sirT 



V r « J 



where 9 n is the n th stored angle of movement direction. 

In the case of flat end milling, all disk elements have the same entry and exit 
angles. However, in ball end milling, each disk element has a different radius as shown 
in FIG. 20. Thus, the entry angle and the exit angle must be calculated in every disk 
element. To calculate the precise entry and exit angles of the k th disk, the difference of 
radius between k m and (k-1) th disks must be larger than the grid size. (If this value is 
smaller than the grid size, a significant calculation error occurs because of an 
insufficient number of nodes.) In this invention, a grouping method which is illustrated 
in FIG. 20 is introduced to solve this problem. The method conducts grouping to make 
each disk group satisfy Formula 44. 

[Formula 44] 



r(ri)-r(n-X)>42g 



where r g (n) is the maximum radius of a disk in n th disk group. The entry and 
exit angles are computed in each disk group, where all the disk elements have the 
same values. 

A flowchart for computing the cutting configuration is illustrated in FIG. 21. 

In FIG. 21, CSS (cutter swept surface) is referenced in an article [Y.C. Chung, 
J.W. Park, H. Shin, B.K. Choi, "Modeling the surface swept by a generalized cutter for 
NC verification," Computer-Aided Design, Vol. 30, No. 8, pp. 587-594, 1998.]. That is, 
CSS refers to a swept surface generated when ball end mill, flat end mill, etc. are 
moved in a straight line from one point to another point in a three-dimensional space by 
an NC code. 

In this calculation process of the cutting configuration, an analytical method 
(CSS calculation method) was used to find the height value (z-value) of cutter swept 
surface in a specific node. The z-value is used to determine whether the target node is 
located in machining region or not (if a height of the target node is lower than the height 
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of cutter swept surface in the target node, the target node is located in the machining 
region.) 

Based on the developed cutting force model, surface error model, and the ME 
Z-map model, an optimum feed rate can be calculated through a feed rate scheduling 
algorithm, which will be described below. 

When cutting force is selected as a constraint variable, it is necessary to 
obtain an optimum feed rate for adjusting a maximum resultant force of the cutting 
force using the cutting force model. Two reference cutting forces as shown in Formula 
45 below are considered using a tool transverse rupture strength presented by tool 
manufacturers. 

[Formula 45] 
RF X = SF • TRS ■ 5 
RF 2 = SF . TRS . S 2 

where RF l represents the reference cutting force considered to avoid breakage of 
the tool shank and RF 2 indicates the reference cutting force to prevent breakage of 
the tool edge. SF means safety factor, which is used to make up for unpredictable 
factors such as cutter geometry error or cutter material variation. TRS means 
transverse rupture strength of the tool material. S { is equivalent cross section area 

and S 2 is cross section area of the edge at the unit disk element. and S 2 are 
calculated as follows: 
[Formula 46] 

S, = RL, • D a 

where D a indicates the height of unit disk. RL^ is the equivalent radius of 
the tool as shown in FIG. 22A; RL 2 indicates rupture length when the tool edge is 
broken in FIG. 22B; a, b, and c indicate rake angle, 1 st flank angle, and 2 nd flank angle, 
respectively; and L, and L 2 indicate the length of the 1 st flank plane and that of the 
2 nd flank plane, respectively. F c is used to indicate the position where tool edge 
breakage occurs. 

RL 2 indicates the distance between E } (x x ,y { ) and E 2 (x 2 ,y 2 ). E x is calculated 



as (F c tan a-F c ) and E 2 is determined as the intersection point between RL 2 and 
the 1 st flank plane or between RL 2 and the 2 nd flank plane. 



2 can 



■r n r i ii_ j et « i i L . A . A tan a(tan a + tan b) „ „ 
If /tt, meets the 1 flank plane, that is, ^ -<F C tana, £ 

(tan"a + l) 

be calculated as follows: 

[Formula 47] 

F r (l + tan 2 a) 

^2 =(77^ — i -tan 6*,) 

(tana + tanZ>) 

ir r»r l xi_ ^st n • i j i « . L x tan a(tan a + tan b) 

If #Z, 7 meets the 2 sl flank plane, that is, — - ^ ->K tana , E, can 

(tan 2 a + l) 

be calculated as follows: 
[Formula 48] 

E-, = ( 5 (£, (tan c - tan Z>) + F c (1 + tan 2 a)) -jc, - tan c + Z, (tan c - tan b)) 

(tan a + tan c) 

Since the reference cutting force is given as the resultant cutting force in the XY 
plane, RF 2 , the normal force on the rake force should be transformed as follows: 
[Formula 49] 



RF 2 '= N d - RF l • ^/{\+B\) 

where B x is the parameter of K f used in Formula 30, and N d means the number of 
disk elements engaged in a workpiece. 

Finally, the minimum value between RF ] and RF 2 is selected as a reference 
cutting force as follows: 
[Formula 50] 

Reference cutting force = Min(RF l , RF 2 ) 

If the reference cutting force is determined using such a method, a feed rate that 
satisfies the reference cutting force can be obtained using a feed rate scheduling 
method based on the cutting force model. 

In the case where surface errors are selected as the constraint variable, the 
optimum feed rate for adjusting a maximum value of surface errors can be obtained 
using the machined surface error model. The feed rate that satisfies the given 
allowable error is calculated from the machined surface error model. 
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[Feed Rate Scheduling Algorithm] 

FIG. 23 is a flow chart of an algorithm used to calculate a feed rate that 
satisfies a reference value in order to adjust a constraint variable (cutting force or 
machined surface error) in the given simulation interval. The cutting force model and 
5 the machined surface error model are used to calculate constraint variable values. The 
algorithm appearing in FIG. 23 can be used for regulation of cutting force and 
machined surface error. 

In the present invention, there is proposed a simulation time reduction method 
for feed rate scheduling. This algorithm first determines a position of a tool path that 

10 generates a maximum value of a constraint variable (cutting force/surface error) during 
one rotation in the given simulation interval, then provides a maximum feed rate at that 
position for adjusting a maximum constraint variable value as a standard value. When 
the constraint variable value is estimated, it takes a significant amount of time to 
calculate tool deflection always from 0 degrees to 360 degrees. To prevent an increase 

15 in the simulation time occurring as the calculation is repeated, feed rate scheduling is 
performed only at a position of a specific tool path. 

The actual uncut chip thickness in the cutting force model and the surface 
error model of the present invention is calculated considering tool deflection and run 
out, and positions passed by previous edges. Although this method increases the 

20 precision in estimating cutting force and surface error, since it uses an algorithm that 
repeats computations to converge on a correct solution when calculating tool deflection, 
the time required for calculation is increased. Accordingly, when estimating a position 
where a maximum value of a constraint variable in a tool path occurs and a position 
angle where a maximum value of a constraint variable in one rotation occurs, a 

25 geometric uncut chip thickness is calculated as shown in Formulae 10 and 32. As 
shown in FIG. 23, if the user sets a specific simulation time interval when performing 
feed rate scheduling, a constraint variable value is estimated along a tool path in the 
given simulation time interval using Formula 10 or Formula 32, which are geometric 
uncut chip thickness model formulae, and after the location at which the maximum 

30 value of the constraint variable occurs is determined, the position angle where the 
maximum value of the constraint occurs is determined at that location. Subsequently, 
simulation is performed at only the position angle where the maximum constraint 
variable value occurs using Formula 12 or Formula 34, to quickly obtain a feed rate that 
satisfied a reference value. 
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The maximum value of the constraint variable has a tendency to increase as a 
feed rate increases in the given axial depth of cut and cutting speed. Using this relation, 
a feed rate corresponding to a reference cutting force can be obtained through a linear 
search. With reference to FIG. 23, when maximum constraint variable values of specific 
5 feed rates and f 2 are CWt and CW 2 , the feed rate of f nex t corresponding to the 
reference value RV, can be approximated using Formula 51 below. 
[Formula 51] 

f = f , {w-cw x ){f 2 -f v ) 

J next Jl CW 2 -CW X 

With reference to FIG. 23, in the case where the constraint variable calculated 
10 from fnext is different from the reference value, the process is repeated whereby 
whichever of fy and f 2 that is closer to f ne xt is designated as the new f lf and f ne xt is 
substituted for f 2 , then the process for obtaining a new f ne xt is repeated. 

That is, following calculation of a constraint variable CW nex t of when the feed 
rate is f ne xt, it is determined using Formula 52 below if the constraint variable value 
15 CWnext when compared to the reference value RV is less than the error limit. When it is 
less than the error limit, the feed rate f ne xt is applied to the NC code. However, when 
this value is not less than the error limit and the reference value RV is greater than the 
constraint variable value CW nex t, the feed rate f 2 is replaced by f nex t. Further, when this 
value is not less than the error limit and the reference value is not greater than the 
20 constraint variable value CWnext, the feed rate U is replaced by f nex ti then the process 
of obtaining f nex t is repeated. 
[Formula 52] 

CW next -RV 

— < Error Limit 

RV 

Since the feed rate is a natural number, an optimum feed rate can be quickly 
25 obtained through such a method. 

Automatic acceleration and deceleration by a controller is realized in 
exponential function form in an A, C interval as shown in FIG. 24. Accordingly, the time 
it takes to reach a designated feed rate after each block starts is slightly greater than a 
designated acceleration and deceleration time T c . In all blocks of an NC code, an 
30 acceleration and deceleration time of approximately T c is required to reach a desired 
feed rate during machining. In manufacturing machines that are being developed in 
recent times, a look-ahead algorithm is applied to somewhat overcome this problem 
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related to acceleration and deceleration. However, in manufacturing machines that do 
not have a look-ahead function, it is necessary to take into account such acceleration 
and deceleration characteristics during machining. Even when a look-ahead function is 
provided, the acceleration and deceleration characteristics cannot be completely 
ignored. 

If a look-ahead function is not supported in a controller, it is not possible to 
obtain a desired feed rate in a block that is smaller than the acceleration and 
deceleration time of the manufacturing machine. Therefore, when the NC code is 
divided to optimize feed rate, it is necessary that the block size be greater than the 
acceleration and deceleration time. In order to perform machining of a block to a 
predetermined feed rate, the size of the block must be at least two to three times 
greater than the acceleration and deceleration time, and the greater the block size the 
longer the predetermined interval of the feed rate becomes. Accordingly, in the present 
invention, although an NC code block is divided when the amount of change in a 
constraint variable value surpasses a reference value, when the size of the divided 
block is smaller than a size determined by considering acceleration and deceleration, 
the NC code block is divided at the position where the block starts to be greater than 
the limited size. 

In the off-line feed rate scheduling method of the present invention described 
above, a feed rate scheduling method is developed based on a cutting force model and 
a surface error model. Therefore, machining precision and productivity can be 
improved during CNC machining using the proposed feed rate scheduling model. 

Although embodiments of the present invention have been described in detail 
hereinabove in connection with certain exemplary embodiments, it should be 
understood that the invention is not limited to the disclosed exemplary embodiments, 
but, on the contrary is intended to cover various modifications and/or equivalent 
arrangements included within the spirit and scope of the present invention, as defined 
in the appended claims. 
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